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The Comment by Amore et al. [gr-qc/0611029| contains a valid criticism of the numerical precision 
of the results reported in a recent paper of ours [Phys. Rev. D 73, 044022 (2006)], as well as fresh 
ideas on how to characterize a quantum cosmological singularity. However, we argue that, contrary 
to what is suggested in the Comment, the quantum cosmological models we studied show hardly 
any sign of singular behavior. 
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The previous Comment [l| presents a criticism of the numerical precision and the physical interpretation of our 
recent results on the quantum dynamics of Friedmann-Robertson- Walker cosmological models with radiation and a 
negative cosmological constant y. 

Concerning the accuracy of the numerical results the criticism isjustified, although this has played no decisive role 
in our analysis, as acknowledged by the authors of the Comment The best numerical results for the energy levels 
that we obtained by means of the Chhajlany-Malnev method are indeed correct up to 20 digits, but the worst are 
correct only up to 5 digits. The lowest precision typically occurs for the intermediate energies in the spectrum. The 
need to deal with a modified potential instead of the exact one, inherent in the Chhajlany-Malnev method, gives rise 
to significant errors, particularly in the case k = —1. Therefore, the numerical method that we used is not uniformly 
precise, contrary to what is suggested in our paper. The "Variational Sine Collocation Method" employed by Amore 
et al. is more uniformly precise and appears to converge more rapidly than the method chosen by us. Thus, the 
powerful numerical tool used in Ref. l| will probably be useful for future research in quantum cosmology. 

As to the physical interpretation, we have some points of disagreement with the authors of the Comment. They 
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criticize us for basing our conclusions in regard to singularity avoidance solely on the behavior of the expectation 
value of the scale factor, and advance some interesting ideas on singularity identification that deserve a closer scrutiny. 
Firstly, they propose a comparison between the classical and quantum probabilities for finding the scale factor near 
a = Q. In the cases fc = or fc = 1 the classical probability of finding the scale factor in the interval (0, e) is given by 

for e small enough, since the potential T4(a) vanishes at a = 0. As to the quantum probability, one has 



r\^{a,t)\'da^\^{e,t)\'e 
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Since Nd and £ are both independent of e, it follows that 

cx |1'(e,i)P ^ as (3) 



Pq{a < e) 



Pci{a < e) 

because the wave functions employed by us vanish at a = 0. Therefore, it is impossible for the quantum probability to 
be greater than the classical probability in a sufficiently small neighborhood of the classical singularity. The authors 
of the Comment claim that this is merely "a consequence of the Dirichlet boundary conditions on the wave function" . 
However, the boundary condition '!'(« = 0, = is a special case of the general boundary condition 

*'(0,t) = a*(0,i) , aeRU{oo}, (4) 

which is necessary to make sure that the Hamiltonian operator be self-adjoint (the prime denotes derivative with 
respect to a). This includes wave functions such that both ^'(a,t) and ^'{a,t) vanish at a = 0. In these cases, as 
e — !■ the quantum probability Pq (a < e) tends to zero even faster as compared to the classical probability Pd (a < e) . 
Thus, there is a large set of states for which the quantum probability of finding the scale factor near a ~ is arbitrarily 
smaller that the corresponding classical probability. Taking Pq{a < e) > Pci{a < e) for all sufficiently small e as a 
singularity indicator, an idea further discussed below, it is clear that there is a large set of nonsingular quantum states. 
Of course, this does not prove that all states are nonsingular, which prevents us from drawing any sharp conclusion. 
The case k = —1 requires some qualifications. The effective potential in this case is 

V;(a) = 12|A|a'^-36a^ (5) 

and a simple sketch of the graph of Ve shows that for positive energies there is no classically forbidden region around 
a = 0. Thus, the statement of Amore et al. jl|] that for fc = — 1 the classical probability vanishes in a sufficiently 
small neighborhood of a = is correct only if the energy £ ~ 12E is negative. One can easily construct wave packets 
with positive mean energy which invalidate the arguments presented in the Comment in the case = — 1 as well. 

The singularity criterion Pq{a < e) > Pci{a < e) suggested by our critics is a welcome attempt to improve our 
means of recognizing a quantum cosmological singularity. Unfortunately, it is not exempt from difficulties. If A; = — 1 
and the energy is negative the classical model is nonsingular (a = cannot be reached) and Pci{a < e) = for e 
sufficiently small. The fact that Pq{a < e) > Pci{a < e) just because Pci{a < e) = should not be taken by itself as a 
symptom that the quantum model is singular, as our critics seem to do. Since classically forbidden regions are quantum 
mechanically allowed, this criterion would force us to conclude that all nonsingular classical cosmological models are 
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singular at the quantum level. This would mean that, instead of contributing to the avoidance of singularities, the 
quantization process creates singularities where classically there was none, which does not sound very sensible to us. 

According to the authors of the preceding Comment, the study of the behavior of the expectation value of the 
scale factor alone is not sufficient to draw sharp conclusions on the presence or absence of singularities. Although 
they may be right, we cannot fail to remark that the expectation value of the scale factor has often been used as a 
singularity indicator in quantum cosmology ^|. Our critics assert that it is more appropriate to study the behavior 
of the probability of finding the scale factor below a given value. Intuitively, if Pq{a < e) ever became close to one for 
very small e the probability density |4'(a, t)\'^ would be sharply peaked near a = 0, and we would have an indication 
of a singularity. But then most likely (a) would be very small too, so that Pq{a < e) appears to be a singularit 
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indicator essentially equivalent to (a). Both our analysis and the more accurate numerical study of Amore et al 
of the dynamical evolution of the models give no hint that such an extreme condition is reached of very small (a) 
or Pq{a < e) close to unity for very small e . This behavior of the wave packets is reflected — not so faithfully as 
we would like, we concede — in the expectation value of the scale factor, and provides evidence that the quantum 
models do not develop singularities. Our critics state correctly that the expectation value of the scale factor can never 
vanish, but this does not mean that "hardly any conclusion can be drawn from the study of (a)" . Our conjecture that 
Pq{a < e) and (a) are essentially equivalent singularity indicators should not be mistaken for a rigorous statement, 
and the proposal of the study of Pq{a < e) as a more reliable singularity indicator deserves further investigation. 

As a matter of fact, it is hard to make clear-cut statements on presence or absence of cosmological singularities in 
the quantum regime because there is no general agreement on what constitutes a quantum singularity. This is one of 
the fundamental problems of quantum cosmology: there is no indisputable definition of a quantum singularity. Owing 
to the lack of a general, physically reasonable and mathematically precise characterization of a quantum singularity, 
all statements on existence or nonexistence of singularities in quantum cosmological models are open to criticism, and 
this is clearly recognized in the conclusions of our paper Q . The suggestion of the authors of the Comment that one 
should look at the probability Pq{a < e) is interesting and worthy of a deeper analysis, although we suspect that it 
also fails to provide an indubitable singularity indicator. 

In short, our critics are right in considering that the evidence we provided was not sufficiently convincing to show 
beyond doubt that the quantum cosmological models we studied are nonsingular. We hope we have supplied more 
satisfactory arguments in this Reply. Regrettably, to our knowledge sharp and conclusive arguments were not available 
when we wrote our paper, and remain so. At any rate, the suggestions made in the Comment seem a valuable effort 
towards a better understanding of quantum cosmological singularities. 
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